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ABSTRACT 

Accretion flows with pressure gradients permit the existence of standing waves which 
may be responsible for observed quasi-periodic oscillations (QPO's) in X-ray binaries. 
We present a comprehensive treatment of the linear modes of a hydrodynamic, non- 
self-gravitating, polytropic slender torus, with arbitrary specific angular momentum 
distribution, orbiting in an arbitrary axisymmetric spacetime with reflection symme- 
try. We discuss the physical nature of the modes, present general analytic expressions 
and illustrations for those which are low order, and show that they can be excited 
in numerical simulations of relativistic tori. The mode oscillation spectrum simplifies 
dramatically for near Keplerian angular momentum distributions, which appear to 
be generic in global simulations of the magnetorotational instability. We discuss our 
results in light of observations of high frequency QPO's, and point out the existence of 
a new pair of modes which can be in an approximate 3:2 ratio for arbitrary black hole 
spins and angular momentum distributions, provided the torus is radiation pressure 
dominated. This mode pair consists of the axisymmetric vertical epicyclic mode and 
the lowest order axisymmetric breathing mode. 

Key words: accretion, accretion discs - black hole physics - relativity - X-rays: 
binaries. 



1 INTRODUCTION 

Observations of quasiperiodic oscillations (QPO's) in the X-ray light curves of black hole and neutron star X-ray binaries 
have motivated the theo ry of relativisti c "d iskoseismolo gy" , the study of hydrodynamic oscillation modes of geometrically 
thin accretion discs fsee IWagonenll999l and lKatoll200ll for reviews). Standing waves with discrete frequencies can exist in 
such discs because the radial profiles of general relativistic test particle oscillation frequencies can create finite regions where 
modes are trapped. Because the disc is considered to be geometrically thin, radial pressure gradients are usually negligible in 
the equilibrium structure of the accretion flow, which is presumed to consist of fluid rotating on nearly circular geodesies. 

However, radial pressure gradients themselves can permit the existence of discrete modes because they produce a disc, 
or torus, of finite extent. At least in the case of black hole X-ray binaries, QPO 's are only observed in the "hard" and 
"steep power law" spectral/variability states (e.g. iMcClintock fc Remillardl bOQfilh states in which the flow is not solely 
composed of a standard geometrically thin accretion disc. It is therefore conceivable that radial pressure gradients may be 
important in trappin g modes, and vario u s gro ups have recently explored the possibility t hat ac c retion tori may explai n 
observed QPO's, e.g. iGiannios fc Spruitl hooj] for low frequency QPO's, iRezzolla etaO te003al) : iKluzniak et alJ \2004) ; 
Lee. Abramowicz. fc Kluzniakl (|2004h for high frequency QPO's. Geometrically thick tori may form in stellar collapse, and 
mode s of dense tori around black holes have also been suggested as a detectable source of gravitational waves fe.g. lZanotti et, all 
l2003h . 

A complete analysis of the spectrum of modes in tori has not yet been done except in the limiti ng case of constant specific 
angular momentum, slender tori in a Newtonian point mass gravitational potential felaesl Il985lh In this paper we greatly 
generalize that early work and present a rather complete analysis of the oscillation modes of non-self-gravitating, polytropic, 
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hydrodynamic tori with small cross-section and arbitrary specific angular momentum distributions. The torus equilibria are 
assumed to be stationary and axisymmetric, and orbit in an arbitrary background spacetime which is itself axisymmetric and 
possesses reflection symmetry. Although limited to slender tori, our work here could be numerically extended to thicker tori 
by tracking all the mode frequencies. Armed with a complete understanding of all the mode frequencies and their physical 
properties, it should be possible to better identify what modes might be responsible for observed QPO's. 

An analysis of torus oscillation modes is also useful as a code check for global numerical simulations of accretion flows. 
In part because of their finite spatial resolution, these simulations generally use tori rather than geometrically thin discs 
as initial conditions fe.g. lDe Villiers. Hawlev fe Krolik1l2003t Icammie. McKinnev fe Totblbooih . It is noteworthy that global 
simulations of non-radiative, m agnetized accretion flow s often produce small, albeit highly variable, hot inner tori near the 
innermost stable circular orbit drlawlev fe Balbusll2002h . 

This paper is organized as follows. In section 2 we derive the basic equations describing the equilibria and linear pertur- 
bations of relativistic slender tori. We then make some general remarks about the types of modes and their classification in 
section 3. In section 4 we present complete solutions for all the modes of relativistic fluid tori in four special cases where this 
can be done analytically. Then in section 5 we show how the low order modes of the general torus can be derived exactly, 
and also present numerical solutions for higher order modes. In section 6 we compare these analytic results with numerical 
simulations of oscillating hydrodynamic tori, and we finally summarize our conclusions in section 7. Readers who are primarily 
interested in applications to QPO's should focus on sections 5 and 7. There we identify a robust new candidate for the modes 
comprising the observed 3:2 ratio of high frequency QPO's in black hole X-ray binaries. We apply our analytic results for 
these mode frequencies to GRO J1655-40 in Figure 9. 

W e remind the reader that slender hydrodynamic tori are violently unstable to the Papaloizou-Pringle instability JPapaloizou fe Pringlj 
Il984l) . Perhaps even more important, in the presence of a weak magnetic field, they are unstable to the development of mag- 
netohydrodynamical turbulence due to the magnetorotational instability ijBalbus fc Hawlevll998l) . Whether or not the hydro- 
dynamic modes we discuss here can survive these instabilities or even exist is a problem which still needs to be investigated, 
and we discuss this briefly in section 7. 



2 SETTING UP THE PROBLEM 

We consider an ideal, hydrodynamic flow in a stationary, axisymmetric spacetime with reflection symmetry. The line element 
may be expressed as 

ds 2 = g tt dt 2 + 2g tr j>dtd(f) + g rr dr 2 + geedd 2 + g^d§ 2 . (1) 

We take the metric to have signaure ( — h ++)■ Because of the assumed symmetries, the metric coefficients g^ v depend on r 
and 9 only, and remain unchanged under the transformation 9 — > 7r/2 — 9. 

We will restrict consideration throughout this paper to flows which are isentropic, and in fact polytropic. Hence the 
internal energy per unit volume e (including rest mass), pressure p, and rest mass density p, all measured in the fluid rest 
frame, are related by 

e — np + p and p = Kp 1+1 ^ n . (2) 
Here n is the polytropic index and K is a constant. We use units where G = c = 1 here and throughout this paper. 



2.1 Equilibrium 

We consider an equilibrium, axisymmetric torus in a state of pure rotation. The only nonzero components of the fluid four- 
velocity are then u* and u^. The angular velocity measured by an observer at infinity is Q, = u^/u 1 , and the normalization of 
the four- velocity u a u a = —1 requires 

A = u* = (—gu - 2£lg t<t , - ti 2 g<i>$) 1/2 . (3) 

Other useful quantities are the specific energy £ and specific angular momentum £ of a given fluid element: 

£ = -u t = {-g tt + 2lg t *-l 2 g^y 1/2 and 1 = -^. (4) 

The spatial r — 9 structure of the torus is given entirely by the relativistic Euler equation, which in this case may be 
written 

^£ 2 (Vg u - 2£Vg t4, + i^g 4 "") = -^-j— Vp. (5) 
Here the operator V may be taken to be either d/dr or d/d9. 
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In the equatorial plane inside the torus is a special radius ro where the pressure has a maximum and fluid orbits on a 
circular geodesic. From equation , the "Keplerian" specific angular momentum £o at this radius satisfies 

(g2-2log* + $g*+) = 0, (6) 

where subscript zero refers to the pressure maximum. 
It is convenient to introduce an effective potential 

U = g tt -2e g t * + %g++. (7) 

A test particle with specific angular momentum £q orbits on a circular geodesic at an extremum of U. This extremum is 
of course located at the pressure maximum ro. Linear perturbations of a test particle about this circular geodesic lead to 
radial and vertical epicyclic oscillations. The squares of the corresponding angular frequencies, as measured by an observer at 
infinity, are given by 



2 El ( 1 d 2 U\ j 2 £$ ( 1 d 2 U 



and "2 > = ^ ^—7^ . (8) 



Al \2g„ dr*J A 2 \2g e 

respectively jAbramowicz et al1l2005() . 

Given our polytropic equation of state ©, the pressure and density can be expressed in terms of an auxiliary function 
/M): 

P = Pof n+1 and p = p f n . (9) 

Consider now a slender torus in which po / po is very small. By introducing local coordinates measured from the equilbrium 
point, 



,1/2 (I 



T 0\ „„J .. - 1/2 ( T/ 2 



and y = g e ' go — , (10) 



ro / \ r 

Abramowicz et al. (2005) expanded equation JSJ to show that 



_ 2 _ g 2 ( g%-e g? de 

Affllio \ g rr dr 



x 2 +wjy 2 }. (11) 



P 2 - SSr ( 12 ) 



Here 

, 2 _ 2(w + l)p 
p A 2 r 2 n 2 

is a dimensionless parameter which determines the thickness of the torus. The slender torus limit corresponds to /3 — > 0. The 
quantities Q r and Qg are the epicyclic frequencies scaled by the angular velocity at the pressure maximum Qq. 

It t urns o ut that the equilibrium function / given by equation 11 U can be written more simply and physically. Following 
ISeeuinl Jl975h , define the spatial vector 

7 = AS 3 [(1 + m) Vg^ - flVg u - iVg^] (13) 

For barotropic con figurations such as the polytropic tori we are considering here, surfaces of constant f2 and £ coincide 
i Abramowiczlll97lh . By taking the curl of equation 0, it is straightforward to show that the vector 7 is then perpendicular 
to these level surfaces, 

7 x VI = 0. (14) 
Even more importantly, the quantity 



A 2 ' A 2 \ g rr dr gee 09 

represents the characteristic squared frequency (as measured by an observer at infinity) of inertial oscillations in the fluid 
due to the presence of an equilibrium specific angular momentum gradient. We show this explicitly later in this paper, but 
the result is not surprising giv en that -y ■ V I is the quantity that appears in relativistic generalizations of the Rayleigh and 
Hdiland local stability criteria l|Segumlll975h . At the pressure maximum, equations @, l|13|l . and l|15^ imply 

k = k S2 = — — 2 ^- . (16) 



£ Al \ g rr dr 
Hence equation li lt may be written as 

/ = 1 - (pr - Ko) - ^ey 2 , (17) 
where x = x/f3 and y = y//3. 
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2.2 Perturbations 

The equations describing the dynamics of linear perturbations of the torus may be derived by linearizing the relativistic 
continuity equation, 

V a (pu a )=0, (18) 
the relativistic Euler equation, 

JT VflU = 2 VnP (19) 

e +p 

the equation of state 10, and the normalization of the four- velocity, about the equilibrium discussed above. Because the 
equilibrium is st ationary and axisymmet ric, we choose all the perturbation variables to vary with t and 4> according to 
exp[z(mc6 — tjt)]. Ilpser fc Lindbloml (Il992l) have shown that the resulting perturbation equations can then be combined into 
a single, second order partial differential equ ation in terms of a singl e scalar perturbation variable. In order to follow the 
notation of recent work on oscillations of tori I Abramowicz et aflboOB^ . we define this variable slightly differently as 

(20) 

pAo 

where a = u> — mQ. We deno te this scalar variable by W , as this was the notation first used in the original Newtonian study 
of global instabilities of tori iPapaloizou fc Pringlelll98i) . 

The pe rturbation equation for to ri could be derived by changing variables and taking the appropriate limit of the general 
equation of llpser fc Lindbloml Jl992h . but we sketch the direct derivation here in order to include some intermediate results 
which will be useful. The direct derivation is greatly aided by noting that equations © and 11311 imply 



AaW 

Vp + V(AaW) = AaVW + Wa 

e + p 



^j + £(A 2 n + g t4, )V£ 



Wm 



£ 



(21) 



Using this and equation I I H , it is straightforward to show that the r and 9 components of the Eulerian four- velocity pertur- 
bation are given by 

5u = ir~ v vw + a*{J-k2) {jr-p) { 7(V ^ ■ vw) + aA£W [ a(A2 " + gt ' ] + ~ 9 ^ ve } • (22) 

The azimuthal component is given by 

Su,/, = -^-^A£{al + mm - m)W + A ^f_ ^ (jT^j V£ ' ^ A °VW + W£ [a(A 2 fl + g^) + m^ff* - g^)] V^j (23) 
These equations are completed by the linearized continuity equation, 

\nll)p W ~ ^ ae + mM ~ m )P^ ~ + TZ^TT? ^9) 1/2 P9 rr Su r ] + ^-^JL[(- g y/* P g ee 6u e ] = 0. (24) 

So far, these equations are exact. In the slender torus limit, r and 8 derivatives of W become dominant in equations l|22|l and 
\\2'M . which then imply 

cj 2 dW . .dW , . £$ff rr fd£\ dW 



8u r ~i— J—— , Su ~ i— -, and 5u$ ~ -* ^g r I 5- ) -j— ■ (25) 

erg — Kg Or 00 erg — Kg V Or / ar 

Substituting into the linearized continuity equation, transforming coordinates from r and 9 to the local, scaled coordinates x 
and y, and taking the slender torus limit j3 — > 0, we finally obtain 

°0 d frndW\ 8 f rndW\ _2-n-l M , n , ofi x 

^-4dE{ f -dk-) + ^{ f -w) +2naof w=0 > (26) 

where ao = ao/iio = — mf2o)/^o- 

This partial differential equation must be solved together with the boundary condition that the Lagrangian perturbation 
in pressure at the unperturbed surface vanish: 

Ap| /=0 = (<5p + e Q, V a p)| /= o = ; (27) 

where is the Lagrangian displacement vector. Using the gauge choice ^ a u a = Jlpser fc Lindblomlll992l) . the r and 9 
components are given by 

I = i-Su. (28) 
Act 

In the slender torus limit, the boundary condition may then be written as 
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r 



g dfdW 1 dfdW 



K 2 dx dx (To dy dy 



+ 2a W 



(29) 



3 GENERAL REMARKS 

Equations 11711 . I26H . and 1291 1 describe a global eigenvalue problem for modes of the relativistic, polytropic slender torus. It 
is important to note that they are, necessarily, identical in mathematical form to the corresponding equations for a polytropic 
slender torus in a Newtonian gravitational potential with axisymmetry and reflection symmetry. All the interesting effects 
boil down to three characteristic frequencies: the radial epicyclic frequency ui r , the vertical epicyclic frequency uie, and the 
characteristic frequency of inertial modes at the pressure maximum «o. 

The slender torus always admits a trivial zero velocity, zero corotating frequency mode, with W being independent of 
x and y, and a 2 = 0. In the constant specific angular momentum case, this is the mode w hich corresponds to the principal 
mode of the Pap aloizou-Pringle instability when the torus has finite thickness feraeslll985Fl . 

As shown bv lAbramowicz et alJ J2005I) . the slender torus (even in the baroclinic case) also always admits incompressible 
modes corresponding to global oscillations of the entire torus at the epicyclic frequencies of the external spacetime: 

W r = x, ol= Col (30) 

and 

We = y, 5o = &$■ (31) 

The fluid velocity for these modes is spatially constant over a torus cross-section. 

Broadly speaking, the remaining modes are of three general types: surface gravity waves, internal inertial waves, and 
acoustic waves. The last two classes are easily revealed in a WKB analysis of equation 1261 . Assuming an x — y dependence 
proportional to exp{i[J k x (x)dx + J k y (y)dy]}, and taking the high wavenumber limit k = (k 2 + k 2 ) 1 ^ 2 — » oo, we obtain the 
local dispersion relation 

at - (k 2 c 2 + ko)(Jo + klc 2 R 2 Q = 0, (32) 



where 

1/2 

(33) 



/3A r n "° \2(n + l)p 

is a scaled local sound speed. This dispersion relation immediately gives two solutions describing the high wavenumber behavior 
for acoustic and inertial waves: 

72-2 

2 — 2 2 2 ^y^O 

ff = k c s and a = (34) 
k A 

respectively. This justifies our claim that Rq is the characteristic frequency of inertial oscillations in the fluid. 

According to the equilibrium equation 11711 . surfaces of constant density and pressure have perfectly elliptical cross-section 
in (x, y) coordinates. The slender torus equilibrium is therefore perfectly reflection-symmetric about x = and y — 0. All the 
slender torus modes can therefore be assigned a definite parity with respect to each of these reflections, and we note these 
parities throughout the rest of the paper. In general, thicker tori remain reflection symmetric about the equatorial plane, and 
so modes will continue to have a definite y— parity. A well-defined x— parity only exists in the slender torus limit, however. 

In addition to the parities, we will also label modes with non-negative integer quantum numbers j and k which are derived 
from the extreme limits of a constant specific angular momentum slender torus in a Newtonian point mass potential and a 
general slender torus with an incompressible equation of state n = 0. These limits are discussed in sections 4.1 and 4.3 below. 
Physically, j > modes have an acoustic character, and have frequencies which diverge in the incompressible limit. Increasing 
j and k corresponds to shorter spatial scales of variation in the mode eigenfunction on a tor us cross-section. 

The tori we are analyzing in this paper are violently unstable iPapaloizou fc Pringleil984h . This instability is not captured 
by our equations, because its growth rate happens to vanish in the slender torus limit. However, for tori of finite thickness , 
the corotation mode becomes the principal mode of the Papaloizou-Pringle instability jBlaeslE985l : iGoldreich et al.lll986l) . 
In addition, the slender torus limit loses the corotation and Lindblad resonances normally associated with nonaxisymmetric 
spiral waves in discs. 



4 COMPLETE SOLUTIONS IN SPECIAL CASES 

There are four special cases where a complete solution to the mode spectrum of the slender torus can be derived analytically. 
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4.1 Constant Specific Angular Momentum in a Newtonian 1/r Potential 

In a Newtonian 1/r potential, the two epicyclic frequencies are both degenerate with the orbital angular velocity: Q r = u>g = 1. 
If in addition the torus has constant specific angular momentum, then Kq = and equation 11711 implies that isobaric surfaces 
have circular cross-section in the variables x and y. By changing variables to (rj, 9) given by 

x — rj cos 8 and y = rjsind, (35) 

Blaes (1985) showed that the partial differential equation 126H becomes separable and solved the complete eigenvalue problem. 
The eigenfunctions and eigenfrequencies may be conveniently written as 

!cos(fc6>) 
or I (36) 
sin(fc(9) J 

and 

a 2 = - (2j 2 + 2jn + 2jk + nk) (37) 
n 

Here Cjt is a normalization constant, j and k are non-negative integers, and Gj is a j-th order Jacobi polynomial jAbramowitz fc Steeunl 
Il972h . 

The zero frequency, j — k = mode is the corotation mode. The j = 0, k = 1 modes correspond to the global radial 
and vertical epicyclic oscillations of the entire torus, which in this case have frequencies that are degenerate with the orbital 
angular velocity. Other modes with j = have frequencies which are independent of the polytropic equation of state, and 
represent surface gravity waves. Modes with k = and j ^ are compressible breathing modes where the torus cross-section 
expands and contracts. All other modes have a mixed acoustic/surface gravity wave character. 



4.2 Keplerian Limit 

A specific angular momentum distribution which is at the opposite extreme to that just considered is a "Keplerian" distribu- 
tion, where every fluid element in the slender torus has the same specific angular momentum as that of a circular geodesic in 
the equatorial plane. In this case the characteristic frequency of inertial mode oscillations Rq is equal to the radial epicyclic fre- 
quency Q r . The eigenvalue problem in this limit becomes very simple even in general axisymmetric spacetimes with reflection 
symmetry. 

In the Keplerian limit, the equilibrium function f = 1 — uj 2 y 2 only varies with height y. One set of modes then has W 
equal to some linear function of x times some function of height y. For them, equation 126H becomes 

d_ ( dW\ 



r iy y r 1W ) ■ 2 " nir u ' "• (38) 

The eigenfunctions are therefore simply Gegenbauer polynomials in y iAbramowitz &: Stegunlll972l) . 

W = CiT~ 1/2) (cooy) or xC^-^iQey), (39) 
where n y is a non-negative integer which is equal to the number of vertical nodes. The corresponding eigenfrequencies are 
n y (n y + In — 1) 



— 2 — 2 
a = Lo g 



2n 



(40) 



The lowest order, n y = mode is the corotation mode, and n y = 1 gives the vertical epicyclic mode. Higher values of n y 
represent standing, ve rtical acoustic waves. (The eigenfrequencies for the n y = 1, 2, and 3 vertical oscillations have also been 
derived bv lKatol2005l within Newtonian theory.) The choice of 1 or a; for the radial dependence indicates that these modes are 
at least doubly degenerate, and in fact we show numerically in section 5 below that there are many degenerate modes whose 
frequencies are given by equation (I4UP at every n y . In addition, we find that all other modes, apart from the trivial corotation 
mode, have frequencies which approach the radial epicyclic frequency. Thus the Keplerian limit of the slender torus consists 
of many degenerate modes at the two epicyclic frequencies and frequencies corresponding to vertical sound waves, plus the 
corotation mode. 



4.3 The Incompressible Torus 

In the limit n — > 0, the torus is incompressible and admits only inertial and surface gravity waves (the latter being the 
j = modes discussed in section 4.1 above). These modes can be explicitly calculated for incompressible tori in the general 
spacetime. In the incompressible limit, equations 1261 and I29t become 
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el 



Kg 8X 



d 2 W d 2 W „ 
+ -^r = 



dy 2 



and 



dfdw i dfdw 

erg — Kg ox aa; <t ay ay 



= 0. 



For modes with ctq > k , it is convenient to change coordinates from x to 

1/2 

/ <7n — Kn ^ 



so that equation 1411 becomes Laplace's equation: 

8PW_ + ^W_ _ Q 

dx 2 dy 2 

The boundary condition 1421 also simplifies: 



1 dfdW 1 dfdW n Tr , 
vo o\ ox o"o ay dy 



0. 



In terms of x> the equilibrium function / now becomes, from equation 1171 . 

/ = i - n 2 x 2 

where 



-2-2 



K = 



-2/-2 

a [u T 



-2\ 



1/2 



It is well-known that Laplace's equation 1441 separates in confocal elliptical coordinates (u, v), defined by 
X = a cosh u cos v 
and 

y = asinhusinw, 

where u is a non-negative real number and ^ v < 2tt. We choose 



1 

TV 



1/2 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



in order that the surface / = of the torus corresponds to one of the elliptical coordinate surfaces. This has a fixed value of 
u given by 



U = tanh 1 ( — \ = - In ( UJe ^—^ ] . 

\CugJ 2 \LUg-1Z) 

The resulting eigenfunctions are then 

= C~£ cosh(fcu) cos(fcu) and = C^T sinh(fcit) sin(fct;) 



(51) 



(52) 



where the C's are again normalization constants and k is again a non-negative integer. The corresponding eigenfrequencies 
are given by 



-2 



klZuJg 



{Qe + ll) k -{ue-llf 
(we + K) k + {Cog - Tl) k 



(53) 



Our derivation of the eigenfunctions (1521 and eigenfrequencies 1531 assumed that a 2 > R 2 . However, equation 1531 
generally admits solutions with <Tq < kq as well. The reader may easily repeat the derivation for this case, and verify that the 
resulting formula for the eigenfrequencies still turns out to be equation 1531 . 

Because 1Z itself depends on <ro, equation lo-'-it represents a polynomial equation for the eigenfrequency which must be 
solved for every k. The universal corotation mode is given by k = 0, while k — 1 gives the two epicyclic modes. For modes 
with k > 1, (To > Rq represent surface gravity waves, while ctq < Kg represent internal inertial waves. For reference, we give 
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Table 1. Lowest Order Modes of the Incompressible Slender Torus. 



x-parity y-parity k ai Eigenfunction 

+ + 1 

+ 1 Gil x 

i <4 y 



2 + ^ ± m + G? e f - 4«^2]l/2 }/2 

2 ^(4^-34)/^+^ -4) l- ^0£vr^_p_^ 



,._2 -2\/-2 -2\-2 -2/4-2 o-2 \ -2 

2/^7,2 _ <»=2n /f ,-,2 , -,2 _ -.2, t _ (4^- Ko )(^-« )^ -^(4^,-3K )j/ 



explicit expressions of the eigenfrequencies and eigenfunctions of some of the lowest order incompressible modes for general 
angular momentum distributions in Table 1. 

The modes become much simpler in the constant specific angular momentum case kq = 0. Then 1Z = u) r , and so 



-2 7 - - 
(Jq = KUJrUle 



(tug + U! r ) 



LOr) 



(uj e + iO r ) k + (ujg — U) r ) k 



(54) 



There are no internal inertial modes in this case, and these are all surface gravity waves. In the limit of a Newtonian 1/r 
potential, u> r = &$ = 1, and equation 1541 gives the j = surface gravity waves of equation 13711 . 

In the Keplerian limit, all the incompressible modes become degenerate with the radial and vertical epicyclic modes. The 
inertial modes all have frequencies approaching the radial epicyclic frequency, while the surface gravity waves approach either 
the radial or the vertical epicyclic frequencies. 



4.4 Neglecting the Vertical Structure 

In an effort to understand the axisymmetric modes seen in numerical studies of tori. lRezzolla et"al"l i2003bll and lMontero et all 



200411 numerically calculated the eigenfunctions and eigenfrequencies of vertically integrated tori. This begs the question - 
what are the slender torus modes if we ignore the vertical direction? 2 The problem is easy to solve in general, and is 
mathematically similar to the Keplerian limit vertical acoustic waves discussed in section 4.2 above. The eigenfunctions are 
again Gegenbauer polynomials: 



w 



-A \ 1 



(55) 



where n x is a non-negative integer equal to the number of radial nodes. The corresponding eigenfrequencies are 



-2 -2 . /-I 
"0 = K o + [U T 



-2\ 



n x (n x + 2n - 1) 



2ri 



(56) 



Note that n x = 1 gives the radial epicyclic mode. IRezzolla et alJ l)2003bh and lMontero et all l 2004h labeled this mode as an 
/-mode, and higher order modes as oi, 02, etc., i.e. o„ x _i. As they point out, this latter set of modes are physically acoustic 
(or p-)modes, at least in this two dimensional situation. However, as we show below in the next section, some of these modes 
are actually related to surface gravity waves in three dimensions, at least in the slender torus limit. This is true in particular of 
the 01 mode, which with the radial epicyclic mode the y propose might be re sponsible for the observed 3 : 2 commensurability 
in high frequency QPO's in black hole X-ray binaries l|Rezzolla et alJl2003al) . 

The ratio in frequencies of the two-dimensional p-modes to the radial epicyclic mode is given by 



(TO , _ 
— = ±ui r 

0J r 



Hi (l 



2n 



- (n x - 1 



, / nx\ «o 
' V + 2n) lj? 



1/2 



(57) 



For constant specific angular momentum and n = 3, the sequence of frequen cy ratios is approximately 1.528, 2, 2.449, in 
good agreement with the (3/2,4/2, 5/2) sequence that IRezzolla et all i2003bl) claimed was an approximate description of their 
lo w order numerical eige nfrequencies. The number 1.528 is also in excellent agreement with the 01/ f ratio found numerically 
bv lMontero et all (|2004|) (see their Figure 6) to be generic for all slender, vertically integrated, relativistic tori with constant 



1 The reader may be surprised that modes with the same eigenfunction can have different eigenfrequencies, as is the case for the fourth 
row of Table 1. However, the relationships between W and the Eulerian perturbations <5p, eq. 121 H . and Su a , eq. 1251 . depend on frequency, 
so the fluid perturba tions a re not the same. 

2 As lRezzolla et alJ J2003bll point out, one must be wary when comparing polytropic, three-dimensional tori with their two-dimensional, 
vertically integrated representations. In general, the latter do not satisfy a polytropic relationship between vertically integrated pressure 
and density. In the slender torus limit, however, there is an exact mapping between polytropic equations of state in the equilibrium 
structure, with ri2D = "an + 1/2 jColdreich et aljl 986jl. Ignoring the v ertical direction is therefore equivalent to adopting the polytropic 
index ri2D- Note that IRezzolla et alJ l2003tJ l and lMontero et al.N2004l) choose n^-Q = 3 for their numerical results. 
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Table 2. Eigenfrequencies of the Lowest Order Modes of the General Polytropic Slender Torus. 
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specific angular momentum. They also found that the frequency ratio dropped for non-constant specific angular momentum 
(see their Figure 8), at least for slender tori, consistent with the fact that the frequency ratios in equation 157H all approach 
unity in the Keplerian limit. We discuss the frequency ratios of the lowest order modes in three dimensional tori in the next 
section. 



5 THE LOWEST ORDER MODES OF THE GENERAL PROBLEM 

As we noted in section 3, the zero corotation frequency and epicyclic modes are always present for the general torus configu- 
ration. One can find additional modes by following the suspicion that a Taylor expansion solution for W about the pressure 
maximum will generally diverge as the surface / = is approached because of the structure of the partial differential equation 
p^ . A solut ion which is finite there, and therefore able to satisfy the surface boundary condition 129H . will therefore probably 
involve a truncation of the Taylor expansion, i.e. involve only finite polynomials in x and y. This is in fact true of all the 
special case modes considered above. 

Guided by the fact that the modes must be even or odd, we try an eigenfunction of the form 

W = axy, (58) 



where a is some constant. This function is odd in both x and y, and clearly satisfies the boundary condition J22J. ft gives in 
fact two modes, with eigenfrequencies determined by the partial differential equation 1261 : 

So = 2 { Q r + &l ± [(«r + u 2 e f - 4kq^ 2 ] 1/2 1 ■ (59) 

These modes are identical to the odd-odd k — 2 modes of the incompressible torus. The positive square root gives a purely 
incompressible surface gravity wave whose poloidal velocity field from equation (12511 is reminiscent of an X-mode polarization 
gravitational wave. The negative square root gives a purely incompressible inertial mode whose poloidal velocity field represents 
a circulation around the pressure maximum. 

We can go further. Try an eigenfunction of the form 

W = a + bx 2 + cf, (60) 

which is even in both x and y. Here a, b, and c are constants, one of which is arbitrary. Once again, the boundary condition is 
trivially satisfied. Substituting this function into equation 12611 . and equating coefficients of 1, x 2 and y 2 , we end up with three 
linear homogeneous equations for a, b, and c. Setting the determinant of coefficients to zero, we obtain the eigenfrequencies 
of the zero corotation frequency mode (for which b=c=0), as well as 



^0 = 2^4 (2n + + O 2 ) -(n+ ± 



(2n + l)(ue — oj'i) + (n + 1)kq I + ^g{^l — «o) 



1/2 



(61) 



The corresponding values of b and c in terms of a are straightforward to calculate, and can be found in Table 3. 

The frequency corresponding to the upper sign in equation I61H blows up in the incompressible n — > limit, revealing 
that it is in general an acoustic mode. It reduces in fact to the acoustic j = 1, k = mode in the limit of a constant specific 
angular momentum torus in a Newtonian 1/r potential, which we discussed in section 4.1. The velocity field in that case is 
a breathing mode - the torus cross-section expands and contracts. In the Keplerian limit of section 4.2, the mode frequency 
becomes that of the n y — 2 vertical acoustic wave. 

The frequency corresponding to the lower sign in equation l|61|l reduces to the even-even k = 2 mode of Table 1 in the 
incompressible, n — > limit. This is a gravity wave, and the velocity field in this limit is reminiscent of a +- mode polarization 
gravitational wave. 

We summarize the frequencies and eigenfunctions of all the lowest order modes of the general slender torus in Tables 2 
and 3. We also show the poloidal velocity fields for these modes in Figure for a torus with Ro/uj r = 0.5. While we are still 
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Table 3. Eigcnfunctions of the Lowest Order Modes of the General Polytropic Slender Torus. 
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Figure 1. Poloidal velocity fields (Su x ,Su y ) of the lowest order, nontrivial slender torus modes. The torus was assumed to be orbiting 
in a Kerr spacetime with a/M = 0.5 and at a Boycr-Lindquist coordinate radius of ro = 10M. The slope of the internal specific angular 
momentum distribution was assumed to give Ro/u) r = 0.5 and the polytropic index was n = 3. 



quite far from the Keplerian limit, the velocity fields are already showing what happens in that case. The X-mode frequency 
becomes degenerate with ujg in that limit, consisting of opposite vertical oscillations on either radial side of the pressure 
maximum. The breathing mode becomes degenerate with the lowest order vertical acoustic wave, as it's velocity field becomes 
largely vertical. The inertial and +-modes consist largely of radial motions in the Keplerian limit, and become degenerate 
with UJ r - 

It is clear that the procedure we used in this section can be extended to calculate even higher order modes. In the special 
case of constant specific angular momentum tori, the dispersion relation remains quadratic for modes of quite high order, and 
we summarize the frequencies and eigenfunctions of these modes in Tables 4 and 5. 

In general, however, higher order modes result in dispersion relations which are polynomials of cubic and higher order, 
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Table 4. Eigcnfrcquencies of the Lowest Order Modes of the Constant Specific Angular Momentum Slender 
Torus. 
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Table 5. Eigenfunctions of the Lowest Order Modes of the Constant Specific Angular Momentum Slender 
Torus. 
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and should be solved numerically. To do this, we substitute the power-series expansion 

w{x,y) = J2 Wii *y- ( 62 ) 

i,l=0 

in equation (12(il . and equate coefficients of x x y l . We find the following set of algebraic equations 
_/°_ 2 [(* + l)(i + 2)Wi +a ,, - (lu 2 - sg)»(i - 

- LD 2 (i + + 2)W i+ 2,i- 2 ] + (/ + 1)0 + 2)W il!+2 
-(a) 2 - + 1)0 + 2)Wi- 2 ,i+2 - S%1{1 ~ 1)W« 

-2 

-2w _ 2 °_ 2 (o; 2 - «g)iWii - 2nLu 2 e Wu + 2na 2 Wu = 0, (63) 

which can be written as a matrix equation A(ao)W — 0, defining the matrix A. We solve equation I|(i3t with basis vectors for 
i = ii + 2(1 — 1) and I = l\ 4- 2(L — 1), where I — 1, J ma x, L = 1, L max . Here ii = for even x-parity and 1 for odd 
^-parity, and similarly for h. Our numerical technique to solve equation l|63|l is to sweep through frequency space and look 
for zeros of the determinant of the matrix A(ao). 

Figure [5] shows eigenfrequencies found by solving equation (163ft compared to the various analytical results. The solutions 
are greatly simplified in the Keplerian limit, as many different branches converge together. The lowest clump approaches the 
radial epicyclic mode at u p . Inertial modes approach uj r from below in the Keplerian limit and go to zero in the constant 
specific angular momentum limit. Surface gravity waves approach ui r from above in the Keplerian limit. The next grouping is 
around the vertical epicyclic mode at ug. Surface gravity waves approach this frequency from above. The next two groupings 
shown are sound waves with 2 and 3 vertical nodes. The frequency increases as more horizontal nodes are added, giving the 
"fan" away from the Keplerian limit. The lowest order sound wave is nearly independent of Ko, and higher order modes are 
more sensitive. 
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Figure 2. Slender torus eigenfrequencies as a function of — Kq)/w^, which ranges from the Keplerian limit (left) to the constant 
specific angular momentum limit (right). Circles show the numerical solution to equation 1261 , assuming four basis functions for x and 
y. Lines show the exact analytic solutions in Table 2. Vertical sound waves in the Keplerian limit are labeled by the number of vertical 
nodes n y = 2,3 (see eauation l4()l . Labels on the right hand side are for modes in a constant specific angular momentum torus (Table 4). 
The black hole was taken to have a/M = 0.9 and the torus centre was at r = lOiW, giving Co I = 0.604 and wj = 0.910. The polytropic 
index was n = 3. The corotation mode (ctq = 0) is not shown, but is also present for all specific angular momentum distributions. 



5.1 Mode Frequency Ratios 

Based on the behavior of its counterpart in vertically integrated tori. iRezzolla et alJ l)2003aF) and lZanotti et, all 1)20051 ) identify 
the "plus-mode" surface gravity wave as ap-mode, which they label as the oi mode. For compressible tori, this characterization 
is not entirely inaccurate. While the radial compressions present in this mode are always accompanied by vertical expansions 
in three dimensions, these vertical expansions can be much smaller than the radial motions if the torus is very compressible, 
as shown in Figure Q 

It is this mode and the radial epicyclic mode that IRezzolla et all ifcOQSah propose might be responsible for the 3:2 
frequency ratio observed in high frequency QPO's of black hole X-ray binaries. Figure|^]depicts this frequency ratio computed 
from equation (161^ for slender, constant specific angular momentum tori of dif ferent polytropic indice s, orbiting Kerr black 
holes with various spin parameters, as a function of the radius of the torus. As IRezzolla et all (l2003al) found numerically in 
various specific cases, the frequency ratio is indeed close to 1.5. Relativity is essential here, as the ratio drops to 2 1//2 in the 
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Figure 3. Ratio of the "plus mode" surface gravity wave = 0, k = 2) frequency to the radial epicyclic mode frequency, for 

a slender constant specific angular momentum (k.q = 0) torus in a Kerr spacetime at different Boyer-Lindquist coordinate radii. Solid 
curves are for a polytropic index n = 3, while dashed curves are for n = 1.5. From top to bottom, each of the three curves at fixed n 
corresponds to spin parameters a/M = 0, 0.5, and 0.998, respectively. The Newtonian limit of the frequency ratio, 2 1 / 2 , is indicated by 
the horizontal dotted line. 



Newtonian limit of large torus radii (see eg. 1371 with j = and k = 2). As shown in Figure |1J the frequency ratio also drops 
for non-constant specific angular momentum, approaching unity in the Keplerian limit Rq = lu^. 

Our general mode analysis reveals that there is another pair of modes whose frequency ratio is not too far from 1.5: 
the "breathing mode" acoustic wave and the vertical epicyclic mode. Figures |5] and HJ show the ratios of these two mode 
frequencies for a variety of equilibrium parameters. For n = 3 tori in particular, which are probably most relevant to radiation 
pressure supported configurations, the frequency ratio is close to 1.5 even in the Keplerian limit. This is in contrast to the 
plus-mode/radial epicyclic mode pair, whose ratio goes to unity in the Keplerian limit. (The reader should compare Figures 
0|and|S| and note the different scalings on the vertical axis.) 

These are but just two examples of possible mode pairs which may be related to the observed 3 : 2 commensurability 
of high frequency QPO's in black hole X-ray binaries. There are presumably others. In particular, we have only considered 
axisymmetric (m = 0) frequency ratios in this section, but of course m/0 would allow an even richer set of possibilities. 



6 COMPARISON WITH NUMERICAL SIMULATION 

A number of authors have performed axisymmetric numerical simulations of hydrodynamic tor i in order to e xplore their 
oscillation modes. Different techniques have been used to excite the oscillations. For example, IZanotti et alJ ll2005l) tried 
initial radial velocity pertu rbations, initial density perturbations , and perturbations based on the ei genfunctions of the modes 
of vertically integrated tori. lLee. Abramowicz. fc Kluzniakl (|20o4 ^ and lRubio-Herrera fc Lee] J2005alfrl have also tried periodic 
external forcing of the inner regions of the torus at a specific frequency. 

Here we use the analytic solutions from section 5 to excite a spectrum of modes in slender hydrodynamic tori and 
then allow them to evolve numerically for a large number of orbital periods. The goals are to verify that we can indeed 
identify the modes at subsequent stages of simulation and lay the foundation for further numerical study of more general 
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Figure 4. Ratio of the "plus mode" surface gravity wave (+, +,j = 0, k = 2) frequency to the radial epicyclic mode frequency, for a 
slender n = 3 torus in an a/M = 0.998 Kerr spacetime at different Boyer-Lindquist coordinate radii and different values of kq. From top 
to bottom, the curves have kq/w t = 0, 0.2, 0.4, 0.6, 0.8, and 0.99. 



l = . 



torus config urations. The numerical evolution was carried out using the Cosmos++ relativistic magnetohydrodynamic code 
described in lAnninos. Fragile, fc Salmonsonl (12005^ . 

The present simulations are restrict ed to two-dimensional axis ymmetry in order to expedite this work and avoid the 
growth of non-axisymmetric instabilities iPapaloizou fc Pringlelll984h . We use a single-level fixed mesh with a resolution of 
128 2 . We use a Schwarzschild spacetime, and the Schwarzschild coordinate grid covers an angular scale 0.488-7T ^ 9 ^ 0.512-7T 
and has radial boundaries at r m in = 9.7M and r max = 10. 3M. The pressure maximum of the unperturbed torus is located at 
ro = 10. 0M, while its inner radius is ri n — 9.9M. The orbital period at the pressure maximum is T ar b = 199M. The polytropic 
index is n — 3 and the polytropic constant is chosen such that the sound speed at the pressure maximum is c s ,o = 0.0011. 
The specific angular momentum distribution is chosen to be a power law, 

where, for Schwarzschild spacetime, A = (— g u / g^^) 1 / 2 ■ 

We initially excite a set of modes by using the eigenfunctions of Tables 3 or 5. The phase of the oscillation modes 
are initial ly chosen such that the perturbations only need to be applied to the velocity fields and not density or pressure. 
Following IZanotti et alJ (120051) , we tracked the modes by looking at the Fourier transform of the norm of the rest-mass 
density ||p|| 2 = 2(=i£ , ii''y anc ^ ma kmg a comparison with the predicted mode frequencies. As we show, this method 
successfully identifies many of the desired modes, but does a poor job on others. It also reveals power spectrum peaks that 
do not correspond to the initial set of modes that were excited. 

Figure [7] shows a portion of the power spectrum of ||p|| 2 for a simulation of a non-constant specific angular momentum 
torus (lo = 3.95 and q — 1.929, so that Ro/Cj r — 0.5). Each of the eigenfrequencies from Table 2 is labeled. Figure |H| shows 
a similar power spectrum for a constant specific angular momentum torus (I = 3.95) with each of the eigenfrequencies from 
Table 4 indicated. While many of the eigenfrequencies are clearly present in the power spectrum, several are not. In particular, 
the vertical epicyclic mode (01+-) is not obvious, even though it was definitely excited by the initial conditions. We suspect 
that this is due to the fact that this mode produces very little in terms of density changes, and so the ||p|| 2 diagnostic is not 
sensitive to its presence. (The fact that the radial epicyclic mode is still easily detected is probably due to the squeezing of the 
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Figure 5. Ratio of the "breathing mode" acoustic wave = 1, k = 0) frequency to the vertical epicyclic mode frequency, for 

a slender constant specific angular momentum (re 2 , = 0) torus in a Kerr spacetime at different Boyer-Lindquist coordinate radii. Solid 
curves are for a polytropic index n = 3, while dashed curves are for n = 1.5. From top to bottom, each of the three curves at fixed n 
corresponds to spin parameters a/M = 0, 0.5, and 0.998, respectively. The Newtonian limits of the frequency ratios from equation 1371 . 
(8/3) 1//2 for n = 3 and (10/3) 1 / 2 for n = 1.5, are shown by the horizontal dotted lines. 



torus by the increasing tidal field of the black hole during the inward phase of the oscillation.) We also tried using \\p — pi\\ 2 
as a diagnostic, where pi is the initial density field of the torus, and found that the vertical epicyclic mode was much stronger 
in the resulting power spectrum, while other modes were missing. More work needs to be done to investigate more sensitive 
numerical diagnostics for these modes. 

The prominent peaks at f/fo — 0.53,1.40 and 1.80 in Figure Q are likely not true oscillation modes. The nonlinear 
interaction of two waves with frequencies /i and /2, if not forbidden by selection rules (parity, etc.), will produce modulation 
at "combination frequencies" |/i ± The three peaks may then be explained as /io++ — /oi+- = 1-53 — 1.00 ~ 0.53, 
0.53 + /o2++ = 0.53 + 0.88 ~ 1.41, and 2/02++ = 2 x 0.88 ~ 1.76, respectively. Alternatively, we find numerically that there is 
a higher order even-even parity mode with frequency 1.78, implying a near resonance with 2/02++- Lastly, we note that there 

is a chance near resonance /01 h + /02++ — /io++- In the slender torus, this interaction is forbidden by the parity selection 

rule for the radial direction. However, this parity rule is strictly valid only for the infinitely slender torus, and may still occur 
in our numerical simulations. 



7 DISCUSSION AND CONCLUSIONS 

We have presented a rather complete discussion of the oscillation modes of non-self-gravitating, relativistic, polytropic tori 
in an arbitrary axisymmetric spacetime with reflection symmetry, in the limit where the tori are very slender. Non- s lender 
toru s oscillation modes can be ex plored numerically, and some work on this has already been done bv lZanotti et alJ J2005I) 
and iRubio-Herrera fc Leel l)2005bf) . It would also be possible to exa mine the beh avior of modes for thicker tori analytically 
using perturbative methods, generalizing the technique employed bv lBlaeJ l)l985h . 

All the slend er tori considered in this pap er are hydrodynamically unstable to the nonaxisymmetric, Papaloizou-Pringle 
instability (PPI) JPapaloizou fc Pringlelll984l They would also devel op magnetohydrodynamical turbulence in the presence 
of a weak magnetic field because of the MRI fealbus fc Hawlevlll99ct) . Almost all simulations examining oscillations of tori, 
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Figure 6. Ratio of the "breathing mode" acoustic wave = 0, k = 2) frequency to the vertical epicyclic mode frequency, for 

a slender n = 3 torus in an a/M = 0.998 Kerr spacetime (dotted curves) and a Schwarzschild spacetime (solid curves) at different 
Boyer-Lindquist coordinate radii and different values of kq. From top to bottom, the curves have kq/lo,. = 0, 0.2, 0.4, 0.6, 0.8, and 0.99. 
The Ko/cJr = 0.99 curves are of course close to the Keplerian limit (eg. !4Ul with n y = 2), shown as the horizontal dashed line. 



including those presented in this paper, have been axisymmetric and neglected magnetic fields, so both of these instabilities 
have not been allowed to develop. A major outstanding issue is therefore whic h, if any , of th e hydrodynamic torus modes 
survive these instabilities. A preliminary investigation of this has been done bv iFragilg teOQot) . but more work needs to be 
done . At least in radially extende d tori, the PPI can be suppressed by dynamical accretion through the torus inner edge telaesl 
Il987t be Villiers fc Hawlevl Eo02 ) . Moreover, the MRI appears to dominate the PPI when magnetic fields are inc luded, and 
generally leads to a near-Keplerian distribution of specific angular momentum in the flow faawlev fc Balbusl BoC^ ) . If modes 
which are essentially hydrodynamic in nature can exist in the presence of MRI turbulence, then the near Keplerian limit is 
therefore probably of most interest. The turbulence itself may then determine the mode amplitudes, but all of these questions 
remain to be investigated. A systematic understanding of all the types of modes which may be present, as we have done in 
this paper, should aid such investigations. 

Figure |21 shows that all mode families converge to the epicyclic frequencies or vertical sound waves in the Keplerian limit. 
If physical discs are indeed nearly Keplerian, this greatly simplifies seismology as only a few distinct frequencies exist. But 
how close are real discs to the Keplerian limit? Equation (8) can be rewritten in terms of the angular momentum gradient for 
Keplerian motion, dlx/dr, by differentiating equation (6) and substituting into equation (8). The result is exactly the same 
as for re in equation (16), with dl/dr, the angular momentum gradient for the fluid torus, replaced by dlx/dr. The difference 
of these two quantities becomes 



2 

K = 



loAl 



g u r - lg*, 



d (1 
d- r {lK 



0: 



(64) 



whi ch manifestly goes to zero in the Kep lerian limit. 

be Villiers. Hawlev fc Krolikl (§003) present plots of Ik and I (averaged against density over spherical shells) against 
radius for MRI driven accretion onto Kerr black holes. In the region of the "inner torus" (see their Figure 3), they find I has a 
similar slope as Ik, but is smaller by ~ 1 — 10%, the difference going to zero at the marginally stable orbit. Using equation l64l 
this implies {ujI-kI)/ujI ~ 0.01 -0.1 in our notation. Hence Figure 2 shows that low order modes should be quite accurately 



Oscillation modes of relativistic slender tori 17 



le-27 




Figure 7. Power spectrum of ||p|| 2 for a torus with kq/u t = 0.5. Each of the eigenfrequencies from Table 2 are identified. The units of 
the vertical axis are arbitrary, while the horizontal axis has been scaled with the reciprocal of the orbital period at the pressure maximum. 



described by the Keplerian limit, while the deviation from the Keplerian limit increases with number of nodes. For emission 
from an optically thick region, flux variations from the photosphere become smaller as more angular nodes are included due 
to averaging of hot and cold spots. Hence low order modes, which are more nearly in the Keplerian limit, are also likely more 
observable. 

In view of this, the lowest order breathing mode/vertical acoustic wave and the vertical epicyclic mode may be an attractive 
candidate for a pair of modes which have a near 3:2 commensurability for radiation pressure-dominated (n = 3) tori, as this 
persists in the Keplerian limit. This is in contrast to the plus mode and radial epicyclic mode pair, whose frequency ratio 
approaches unity in the Keplerian limit. Identifying the vertical epicyclic mode with the lower of the two observed frequencies 
also allows the black holes to have any spin parameter, in contrast to the case where the lower of the two frequencies is 
identified as t he radial epi cyclic mode. This is illustrated in Figure [5] for the case of GRO J1655-40. This source has a mass 
of ~ 6.3 M p, lO rosz 20 031) . and models which identify the 300 Hz QPO as the axisymmetric radial epicyclic mode require 
a/M > 0.9 iRezzolla etld1l2003al: lTorokll2005ll . The same is true for diskoseismology models which identify the 300 Hz QPO 
with an axisymmetric "q-mode" , which has a frequen cy below the radial epicyclic frequency within the mode trapping region 
fe.g. IWaeoner. Silbergleit fc Ort ega-Rodriguezll2 00lt) . This is in conflict with recent continuum spectral fitting results which 
suggest a lower black hole spin CS hafee et alJl2005T) . If the 300 Hz QPO is instead the axisymmetric vertical epicyclic mode, 
any spin is possible. Many other mode pairs would probably also work, especially if we consider nonaxisymmetric (m 7^ 0) 
modes. However, one still needs to explain why the QPO's occur at the actual observed frequencies, i.e. what singles out a 
preferred radius for the torus in this source. 

More work needs to be done to investigate how torus oscillation modes might manifest themselves in X-ray light curves 
and power spectra. Some work has been done based on tracking photon geodesies from the torus to the observer l)Bursa et alJ 
120041 : ISchnittman fc RezzollalboOflf) . but additional work needs to be done on the photon emission mechanism of the tori 
themselves. The difficulties that some diagnostics have in identifying torus oscillation modes in numerical simulations, even 
when they are present in the initial conditions as discussed in section 6 above, should serve as a warning to numericists 
searching for QPO's in their simulations. The X-ray observations generally have nothing to do with these diag nostics. 

I t remains unclear whether tori are appropriate models for the flow ge ometry for the "steep po wer law state" iMcClintock fc Remillardl 
I2OO5I) where high frequency QPO's are seen in black hole X-ray binaries. iKubota fc Done! i20o4) suggest that the steep power 
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Figure 8. Similar to Figure |7| but for a constant specific angular momentum torus with each of the eigenfrequencies from Table 4 
identified. 



law state consists of a standard geometrically thin disc sandwiched between a magnetized corona which has sucked up much of 
the accretion power out of the disc. Oscillation modes that are analogous to the torus modes with vertical reflection symmetry 
might exist in such a magnetized corona, because the vertical velocities vanish in the equatorial plane, where the high inertia 
disc is located. We hope that our comprehensive hydrodynamic mode analysis may act as a foundation for investigation of 
these and the other questions noted above. 
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